Powder Technology 269 (2015) 566-576 



ELSEVIER 


Contents lists available at ScienceDirect 

Powder Technology 

journal homepage: www.elsevier.com/locate/powtec 



Effect of lifter shape and operating parameters on the flow of materials in 
a pilot rotary kiln: Part II. Experimental hold-up and mean residence 
time modeling 

A.S. Bongo Njeng a b , S. Vitu a , M. Clausse a,c , J.-L. Dirion b , M. Debacq a * 

a Conservatoire National desArts et Metiers, Laboratoire de Genie des Precedes pour i'Environnement, i'Energie et la Sante (EA21), 2 Rue Conte, 75003 Paris, France 
b Universite de Toulouse, Mines Albi, CNRS UMR 5 302, Centre RAPSODEE, Campus jarlard, F-81013 Albi Cedex 09, France 
c ESIEE Paris, 2 Boulevard Blaise Pascal Cite DESCARTES, BP 99 93162, Noisy le Grand Cedex, France 


ARTICLE 


INFO 


A B S T R 


C T 


Article history: 

Received 20 January 2014 
Received in revised form 30 April 2014 
Accepted 26 May 2014 
Available online 4 July 2014 


Keywords: 
Rotary kiln 
RTD 
MRT 
Hold-up 

Lifting flights 


Experiments were carried out on a pilot scale rotary kiln equipped with lifters at room temperature to investigate 
the effects of the kiln slope, rotational speed, mass flow rate of materials, and exit dam height on the hold-up and 
the mean residence time (MRT). The MRT was determined from the residence time distribution measurements 
as detailed in Part I of this work. Two granular solids having different properties were used: sand and broken rice. 
Furthermore, two shapes of lifters were compared to determine the influence of lifter geometry: straight lifters 
(SL) and rectangular lifters (RL). A new model to predict the MRT was established by means of a dimensional 
analysis. The correlation not only gave good agreement with the experimental data from the present study, 
but also demonstrated good predictive performances when applied to published experimental data of other 
kilns; the model is applicable for inclined kilns that process materials in cascading (tumbling) motion, whether 
or not equipped with lifters or fitted with dams at the outlet end. 

© 2014 Elsevier B.V. All rights reserved. 


1. Introduction 

Rotary kilns are gas-solid reactors widely used in mineral process 
applications and other processes applied to specific granular materials. 
Due to their extensive use in industry, there have been several studies 
aimed at modeling the transport of the solids through the kiln cylinder 
[1-13], which is usually equipped with lifters [14-17], Most of these 
studies attempted to understand and predict key parameters such as 
the hold-up, the mean residence time (MRT) or, as shown in Part I of 
this study [18], the residence time distribution of solid particles. 

Most of the substantial earlier body of scientific literature that exists 
in the field assumed no axial dispersion through the kiln, such that the 
time of passage r was analyzed instead of the mean residence time t. 
Others simply conflated the two and assumed the mean residence 
time to be equal to the time of passage. Indeed the time of passage is 
very easy to determine; it results from the ratio of the weight of the 
kiln hold-up to the mass flow rate. However, it should not be forgotten 
that the time of passage might fail to reflect the flow of solids, especially 
when there is some axial dispersion, which is increased when using 
lifters. 


One of the earliest equations was developed by Sullivan et al. [ 1 ] for 
the calculation of the time of passage of particles in rotary kilns which 
do not have lifters, as follows: 


_ unVe 

SDN 


x factor 


(1) 


where 6 is the angle of repose, S is the kiln slope, N is the rotational 
speed, L and D respectively are the length and internal diameter of 
the kiln, and factor is a parameter accounting for the operating condi¬ 
tions and is equal to one for a simple kiln without obstructions or 
constrictions. 

Chatterjee et al. [5] used a semi-empirical correlation to predict the 
mean residence time of solid particles. Their model was built on a di¬ 
mensional analysis taking into account parameters influencing the 
flow of materials. However, some parameters were missing such as 
the exit dam height, which did not appear in the list of parameters, 
though its effect was studied in their experimental matrix. By applying 
conditions of dimensional homogeneity, they deduced the following 
correlation: 
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where F is the volumetric feed rate, k is a constant equal to 0.1026 and 
exponents [a, e, c} are evaluated as {1.054, — 0.981, —1.1} [6]. It can 
be noted in this correlation that the effect of the kiln slope and the 
angle of repose of solid particles are strongly linked. 

In the literature, models for the mean residence time or time of 
passage, developed for kilns both with and without lifters, range from 
relatively simple empirical equations to fully mechanistic models. How¬ 
ever, we will focus here on some of the simple empirical models devel¬ 
oped for kilns equipped with lifters, since the so-called mechanistic 
models are not suitable for control purposes [4], 

Among the earliest models developed for rotary kilns equipped with 
lifters is the relation of Prutton et al. [2] to predict the time of passage: 


(3) 


where V s is the volume of solids inside a lifter, k is a constant depending 
on the number and shape of lifters, and m is a constant depending on 
the gas flow direction and the properties of the materials. This model 
does not consider the effect of any changes in the flow rate of the 
solid particles and/or gas through the rotary kiln. Moreover the con¬ 
stants k and m must be determined by fitting the model to experimental 
data from the rotary kiln in use, thus making the model impossible to 
use for design purposes. 

Certainly the most commonly used model, is the correlation of 
Friedman and Marshall [19], which is derived from the model of 
Sullivan et al. [1 ] with an additional term to allow for air drag on the 
solid particles. The Friedman and Marshall [19] equation was obtained 
through the study of hold-up of a range of materials under overloaded 
conditions: 

, [0.3344 0.6085C1 


where d p is the average particle diameter, and F and G are respectively 
the feed rate of the solid particles and the flow rate of the gas. The 
second term in this equation expresses the air drag; thus the negative 
sign is used for counter-current flow and the positive sign is used for 
co-current flow. Compared to the model of Sullivan et al. [1], this 
model may account for variations in solid and gas feed rates. However, 
by increasing the solid and gas feed rate by proportional amounts, the 
predicted residence time will remain constant. Furthermore, the geo¬ 
metric features of lifters are not taken into account in this formulation. 

Shahhosseini et al. [4] modified the model of Friedman and Marshall 
[19] in order to better represent the dynamics of the system, the main 
objective being to determine the hold-up through the use of the reten¬ 
tion time of sugar in the rotary kiln. An additional term was introduced 
to account for changes in the flow rate of the solid particles under zero 
gas flow rate conditions: 


(5) 


This model attempted to address some of the shortcomings of the 
Friedman and Marshall [19] model, but still did not consider internal 
fixtures of the kiln. 

Alvarez and Shene [20] presented an empirical relationship for 
residence time estimation, ft was derived from experimental data 
concerning solid particles of biological and mineral origin. Their correla¬ 
tion requires up to six constants and five other exponents, leading to a 
total of 11 parameters to be determined: 
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where p butk is the density of the bulk materials. This formulation allows 
prediction of residence times even when the kiln is horizontal and may 
apply to a wide range of solid materials. However, due to the large 
number of model parameters there might be a need to fit the model 
with the system experimental data prior to its utilization. 

Recently Thibault et al. [21 ] proposed two empirical models for 
the prediction of the mean residence time with the use of a data¬ 
set of different types of solids. The models consist of three functions 
as follows: 


t = p p f,(\/FN)f 2 (S)f 3 (G). 


(7) 


The three functions/i,/ 2 and/ 3 account for the influence of the prod¬ 
uct FN, the slope of the rotary kiln, and the impact of the gas flow rate on 
the mean residence time. The definition of f 2 and / 3 as first order func¬ 
tions in the models makes it possible to calculate the MRT for horizontal 
kilns with no gas flow rate. The main difference between the two 
models proposed is the definition of the function/) which can be either 
a second-order polynomial with the inverse of FN or a simple function 
with F and N each with an exponent. 

This paper is concerned with the study of the effects of usual operat¬ 
ing parameters, shape of lifters, and type of materials on the mean 
residence time of solid particles as well as on the kiln hold-up. A semi- 
empirical model obtained from dimensional consideration of the 
mean residence time of solid particles is also presented. Compared to 
the other published correlations, this new one should be applicable to 
a wide range of kilns and operating conditions, and could also be useful 
for design purposes. Therefore this paper will compare the predictions 
obtained with this model when applied to experimental data from the 
literature to the corresponding experimental MRT, and will also present 
the predictions of published models using the present experimental 
data. 


2. Materials and methods 

2.1. Apparatus and materials 

The apparatus and materials have been described in Part I of this 
work. Therefore only their main characteristics are given here. The rota¬ 
ry kiln is a 0.101 m diameter and 1.95 m long tube, which can be tilted to 
a maximum angle of 5°. Its rotational speed can be adjusted by a driven 
motor between 2 and 12 rpm. The flow rate of the solid particles was ad¬ 
justed by regulating the rotational speed of a screw feeder. Different 
dam heights were used as given in Table 1. In order to study the effect 
of internal flights, 4 straight one-section lifters of 10 mm (inside length), 
referred to as straight lifters (SL), and 4 two-section lifters of 10 mm 
(inside length), referred to as rectangular lifters (RL) due to their right 
angle cross section as shown in Fig. 1, were used. 

Quartz sand (average size 0.55 mm) with a narrow distributed size 
fraction and broken rice (average size 3.8 mm in length, 1.9 mm in 
diameter) were employed as experimental materials. The physical 
properties of these materials are listed in Table 2. 


Table 1 

Details of the circular dams at the kiln outlet end. 


Exit dam height Open diameter Diameter covered Area covered 

[mm] [mm] [%] [%] 

0 101.3 0.00 0.00 

13.5 74.3 26.67 46.20 

23.5 54.3 46.40 71.27 

33.5 34.3 66.14 88.54 


SS+1 
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22. Experimental technique 

The experimental mean residence time was determined from the 
residence time distribution (RTD) measurement data. The MRT is de¬ 
fined by [22]: 

f fr-gWi>M 

J C(t)dt ^^ 5 C(ti)Mi 

where C(t,) represents the tracer concentration at the kiln exit end at 
the discrete time t,-, i = {1,2,3.N s }, and Aq the associated sampling 

The kiln hold-up was obtained after the RTD measurement. The 
feeding system was disabled, and then the solid materials within the 
kiln were discharged into the recovery tank and weighed. Thus, it was 
possible to assess the time of passage of particles in the rotary kiln by 
dividing the weight of solid materials measured by the mean flow rate 
at steady state. 

Part I of this work provides further details on the whole experimen¬ 
tal procedure. 

3. Results and discussion 

Table 3 gives the results of the hold-up and the mean residence time 
for a number of experiments where sand and broken rice were operat¬ 
ed, with either straight, rectangular or no lifters. The details of these 
runs are also given. It must be pointed out that the average mass flow 
rate measured was within a margin of ±0.05 kg-h _1 of the reported 
values. Note that in runs where the kiln was equipped with lifters, 
when the steady state was achieved, the kiln was either design-loaded 
or over-loaded, i.e. there were always enough solids at the bottom of 
the kiln to fill lifters to their maximum capacity. Repeated experiments 
are not reported in Table 3. This part of the paper will discuss the effects 
of kiln rotational speed, lain slope, mass flow rate, exit dam height, lifter 
shape and bulk material properties on the fractional volumetric hold-up 


Table 2 

Physical properties of materials. 

Materials Shape p Mt Propped d p 0 

[kgjn~ 3 ] [kg.m- 3 ] [mm]|f 
Sand Nodular 1422 1543 0.55 39 

Broken rice Cylindrical 889 934 3.8 x 1.9 36 


and the mean residence time. Next, a model for the prediction of the 
MRT will be proposed. 


3.1. Hold-up, volume fraction 

As shown in Table 2, the particulate solids used as bulk materials 
have varying properties. It was convenient to use the fractional volu¬ 
metric hold-up (or elsewhere the degree of filling) for the analysis of 
the results. Fractional volumetric hold-up is defined as the ratio of the 
volume of the solid loading to the volume of the kiln cylinder. The 
volume of solids was determined by dividing the hold-up of the solid 
particles by the bulk density. The effects of operating parameters on 
the volumetric hold-up are illustrated in Fig. 2. Note that only experi¬ 
ments with lifters are presented in Fig. 2. At given operational condi¬ 
tions, experiments processing broken rice filled a higher volume than 
those with sand, even though the values of the hold-up weight were 
very close. This is a direct consequence of the materials' differing bulk 
density properties. It is also important to note the effect of lifters: the 
experimental results show that the presence of lifters significantly 
increased the hold-up (see Table 2). The results also show that the use 
of rectangular lifters slightly increased the hold-up within the lain com¬ 
pared to the use of straight lifters. This may be related to the lifters' 
hold-up capacity, which is higher in the case of rectangular lifters. 
Hence the loading of a kiln will be partly dependent on the presence, 
size and shape of the lifters as previously suggested by Prutton et al. 
[2] and Afacan and Masliyah [23]. 

The classical operational parameters also significantly influence the 
volumetric hold-up. The fractional volumetric hold-up sharply de¬ 
creased at smaller rotational speed values before slowing down at 
higher speed values, as if a bottom had been reached (see Fig. 2a). In¬ 
deed the hold-up must be at the lowest possible level when the motion 
is starting to centrifuge. This usually happens at very high speed. The 
fractional volumetric hold-up gradually fell as the kiln slope increased 
(see Fig. 2b). On the contrary, the volumetric hold-up steadily increased 
with the feed rate (see Fig. 2c). This is not obvious but still can be under¬ 
stood if one extrapolates the work of Saeman [3] followed by Descoins 
et al. [9] on the bed depth profiles. These studies show that in a kiln 
without lifters any increase in the feed rate will result in an increase 
in the bed depth and hence an increase in the fractional volumetric 
hold-up. Even if their models are not directly applicable to kilns 
equipped with lifters, the results remain the same with respect to the ef¬ 
fect of the flow rate on the volumetric hold-up. It was observed that the 
volumetric hold-up increased slightly with increasing exit dam height 
(see Fig. 2d). The higher the dam height, the higher the kilning bed 
depth, which explains the obtained results. 
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Table 3 

Experimental results: mean residence time and hold-up. 

Operating conditions RL SL NL 


N S Jtf h MRT HU MRT HU MRT HU MRT HU MRT HU 

[rpm] [°] [kg.h '] [mm] [min] [kg] [min] [kg] [min] [kg] [min] [kg] [min] [kg] 


2 2.5 

2 2.5 

2 2.5 

2 2.5 

2 2.5 

2 2.5 

2 2.5 

2.5 2.5 

3 2.5 

4 2.5 

5 2.5 

2 0.68 

2 1.3 

2 1.9 

2 2.5 

2 2.5 

2 2.5 


23.5 

23.5 

23.5 

23.5 

23.5 
235 
235 
235 
235 
235 
235 
235 
235 
235 

0 

13.5 

33.5 


43.4 1.953 48.6 

67.4 2.953 65.5 

32 3 1.391 39.7 

25.1 1.033 30.0 

20.6 0.849 25.5 

18.2 0.704 22.2 

16.2 0.653 20.5 

31.8 1.305 41.9 

26.4 1.137 35.5 

20 0.858 24.9 

14.6 0.625 18.9 

56.9 0.573 63.9 

483 1.028 54.8 

44.5 1.415 52.1 

39.6 1.685 41.3 

39.1 1.713 42.9 

51.1 2.168 59.6 


1.959 39.2 

2.808 57.9 

1.536 

1.143 21 

0.973 

0.831 15.2 

0.774 

1.631 

1.333 24.9 

0.937 18.9 

0.734 14.5 

0.676 54.9 

1.105 41.6 

1.541 42.2 

1.646 33.5 

1.720 

2.489 48.6 


1.736 44.8 

2.505 62.7 

35.8 

0.932 26.7 

22.2 

0.684 19.3 

17.1 

39.8 

0.919 30.6 

0.705 23.2 

0.625 16.6 

0.570 59.7 

0.967 51.7 

1.276 49.1 

1.478 37.0 

2.122 57.5 


1.886 29.5 


1.508 

1.124 18.0 

0.892 

0.791 13.4 

0.736 


0.918 

0.693 

0.670 

1.085 

1.503 

1.539 

1.619 

2.478 


34.9 


42.1 


1.214 

0.724 

0.532 

0.774 

0.758 

1.742 


3.1.1. Time of passage 

The time of passage can be calculated from the hold-up measure¬ 
ments. as follows: 



Most of the previous studies focus on this simple relation to investi¬ 
gate the transport of solid particles through a rotary kiln, sometimes 
conflating it with the MRT. as previously stated. In the present study, 
therefore, we compared the time of passage and the MRT results over 
the whole experimental campaign as shown in Fig. 3. It must be said 
that, not surprisingly, quite good agreement was found between the 


(a) Variation of hold up volume 
fraction with kiln rotational speed 



▼ Sand-RL 

t 

# Sand-SL 
A Rice-RL 

O Rice-SL 

A 

• ° A 

6 A 

* * 
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o 1 -*-*-*-*- 1 -*- 

0 2 4 6 8 10 12 14 

Kiln rotational speed [rpm] 


(b) Variation of hold up volume 
fraction with kiln slope 


T Sand-RL 
• Sand-SL 
A Rice-RL 
O Rioe-SL 


* * ft 

3 4 

Kiln slope fj 


(c) Variation of hold up volume 
fraction with mass flow rate 


(d) Variation of hold up volume 
fraction with kiln exit dam height 



Feed rate [kg.h -1 ] Exit dam height [mm] 


Fig. 2. Influence of operating parameters (N, S, Mand h) on the hold-up volume fraction for the flow of sand and broken rice, when the kiln is equipped with either straight or rectangular 
lifters (including results from repeated experiments). 
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i 



Time of passage [min] 


Fig. 3. Comparison of experimental MRT and time of passage of solid particles. Solid lines 
are ±20% margins. 


MRT and corresponding time of passage in general. This is because 
the dispersion coefficients observed in the system are quite low, as 
shown in Part I of this work. Contrary to the MRT, this method does 
not offer insight into dispersion related issues, which are linked to the 
residence time distribution of the solid particles. This is why some 
deviations are still observable. The time of passage estimates the mean 


(a) Variation of MRT with 
kiln rotational speed 


4 

o 
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(c) Variation of MRT with 
mass flow rate 
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Fig. 4. Influence of operating parameters (N, S, Mand h) on the MRT for the flow of sand and 
results from repeated experiments). 


residence time only with a certain degree of accuracy depending on 
the installation. 


32. Mean residence time 


32.1. Influence of operating conditions 

The mean residence times of sand particles were found to be higher 
than those of broken rice particles, mainly due to the effect of gravity 
and their differing angle of repose. For a given type of material and 
fixed operating parameters, the MRT obtained while using rectangular 
lifters was higher than that obtained while using straight lifters, due 
to the differing loading capacity of the lifters. These trends, as shown 
in Fig. 4, were previously discussed in detail by Bongo Njeng et al. 
[18], Note that only experiments with lifters are presented in Fig. 4. 

The MRT significantly decreased by about 69% as the rotational 
speed increased 6 times, as shown in Fig. 4a. These results also agreed 
well with the findings of Abouzeid et al. [7] while varying the rotational 
speed over a wide range of values. Abouzeid et al. [7] found in their 
experiments that the MRT declined when increasing the rotational 
speed up to 20 rpm. With a further increase in speed the MRT remained 
approximately unchanged between 40 and 80 rpm. Above 80 rpm, 
when particles were cataracting, they found that the MRT significantly 
increased. Fig. 4b and c show that the MRT decreased with the kiln 
slope and the mass flow rate. But it must be further specified in the 
case of the inclination that there was a sharp decrease between 2 and 
4° followed by a slow decrease for higher slopes, whereas the MRT 
slowly decreased when the mass flow rate increased from 0.68 to 
2.5 kg-h _1 (see Fig. 4c). This tendency is similar to the one observed 
by some other researchers [5-7,24], When increasing the exit dam 
height the MRT slightly increased (see Fig. 4d) as observed previously 


(b) Variation of MRT with 
kiln slope 
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Table 4 

Calculated values of exponents and constant k for the MRT ( 


Exp, data_k_ a _ ft _ y _ 8 _e_ s _ £ _ rj\6]_ 

Sand 0.1363 0.0508 - 0.4008 0.8749 - 0.9814 0.8115 - 4.5285 0.7723 1.1 

Rice 0.0792 -0.0218 -03387 0.8749 -12277 0.8184 -8.0175 0.7723 1.1 

Sand & rice 0.2611 0.0842 - 03649 0.8749 -1.1243 0.8350 - 5.5283 0.7723 1.1 


by Chatterjee et al. [5], However, in our experiments, a maximum in¬ 
crease of only 5% was observed on the MRT by using a 13.5 mm exit 
dam height compared to the case where no exit dam was fitted at the 
kiln outlet end. It has been observed that these two experimental condi¬ 
tions, all in all, led to similar bed profiles. The only difference which can 
be mentioned is that in the first case the bed depth was kept constant till 
the outlet end with the use of the dam, while in the second case it 
started displaying a slope assumed to be close to the angle of repose ap¬ 
proximately 20 mm from the exit. Hence the design of a kiln and its exit 
dam should give careful consideration to the aspects mentioned here to 
optimize the production objectives. For a dam to have a significant 
impact on the transport characteristics of the materials, its height 
must be higher than the bed depth observed when the kiln has a clear 


32.2. Model 

The main purpose of this model is to come up with a prediction of 
the MRT of solid particles defined as a function of the kiln design, the 
properties of the bulk materials, and the operating parameters, namely 
the kiln rotational speed, the kiln slope, the mass flow rate and the exit 


dam height. For the purpose of application to other kiln design and 
operating conditions, and thus to meet scaling-up issues, the model 
for the prediction of the MRT is based on a dimensional analysis as 
described below. 

322.1. Dimensional analysis. In dimensional analysis, a certain function¬ 
al form is assumed to relate the varying parameters. All the possible 
parameters are commonly present with their own exponent. In this 
study the main operating, geometrical, and physical parameters which 
can affect the MRT are listed below along with their dimensions: 


N Rotational speed T -1 

M Mass flow rate MT~' 

S Kiln slope 

D open Effective exit open diameter L 

Sujt Cross section of materials in lifters L 2 

L Length of the kiln L 

D Internal diameter of the kiln L 

p bu ik Bulk density of materials ML -3 

p tapped Tapped bulk density of materials ML -3 

9 Angle of repose of the bulk materials 

g Gravitational acceleration LT 2 


(a) Experiments with sand 








Fig. 5. Comparison of the experimental MRT while using a) a bulk of sand and b) a bulk of broken rice with predicted values from Eq. (13) using the 3 sets of parameters as giv 
Solid lines are ±20% margins. 


l Table 4. 
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The analysis involves the fundamental dimensions MLT: mass, 
length and time. It follows from the above that the relation being sought 
can be summarized as: 


x. AfrLP 2 (n 2 D) 



F(t, N, M, S, D open , Sjjft, I, D, Pb^, 6, g). (10) 

It is convenient to rewrite Eq. (10) so that the dependent variable f 
appears directly. Whatever the function F may be, Eq. (10) can be 
written in the form: 

F(N, M, S, D open , S lift , L, D, p bulk , p tapped , 6, g) ■ i = 1. (11) 


Assuming that Eq. (11) is dimensionally homogeneous, the parame¬ 
ters of function F must be subject to the power dependence. According 
to Buckingham's theorem, the formula connecting the set of parameters 
can be expressed as a function of dimensionless arguments and each of 
these arguments is a dimensionless product of a combination of the 
original parameters raised to integer power. Having in Eq. (11) a num¬ 
ber of variables r = 12 and a number of fundamental units n = 3, the 
number of dimensionless groups p is, according to Buckingham's second 
theorem: 

p=r-n= 12-3 = 9. 


A collection of 9 dimensionless groupings formed from the parame¬ 
ters in Eq. (11) were defined as follows: 


tM JN 2 P Dopen es ft 4V p bulk il 

Pbu&^-D 2 [ £ D p bulk NLD 2 jrD 2 Ptapped Dj 


(12) 


In the literature, some of these dimensionless expressions have been 
defined, i.e. the Froude number, the length-to-diameter ratio, or the (in¬ 
verse of the) Hausner ratio.is the fractional open diameter at the kiln 
outlet end. Sf estimates the fraction of tube cross section not covered by 
solid particles in loading lifters, where S (1 j t = "if—" , ‘*r l S fl0 rfi/t with nujt 
being the total number of lifters and S horHjt being the section of materials 
in a lifter at horizontal position. When there are no lifters fitted inside 
the kiln tube, S hor/ ,/ t = 0 and ^ = 1. The two dimensionless groupings 
an d p % r) i are very similar to those found by Chatterjee et al. [5] 
in their correlation. However, in this work, contrary to what Chatterjee 
et al. [5] proposed, the bulk materials' angle of repose and the slope of 
the kiln are taken as separate factors rather than combined as an angular 
ratio (see Eq. (2)). Eq. (12) can be written as specified by Buckingham's 
theorem so that an expression of the MRT is obtained from the defined 
dimensionless numbers as given in Eq. (13). 

The variables in Eq. (13) are expressed in the international 
system of units, therefore the mean residence time obtained is in 
seconds. The exponents a, /3,.... g and the constant k are to be deter¬ 
mined with the use of experimental data. The length-to-diameter 
ratio was not varied in our experimental matrix, therefore the expo¬ 
nent r; was fixed to the value of 1.1, which was determined by Chat¬ 
terjee et al. [6] when varying the dimensionless group t from 3.3 to 
16. To determine the other exponents and the constant k, a MatLab® 
script was written. The program sought to determine the parameters 
by minimizing the sum of square errors between the experimental 
and predicted values of the MRT. The program was first run with the 





Fig. 7. Comparison of experimental MRT from Colin 


il. [28] with predicted values from Eq. ( 13) using the 3 sets of parameters as given in Table 4. ! 
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(a) Sai et al. [27] 



(b) Chatterjee et al. [5] 





Fig. 8. Comparison of experimental MRT from a) Sai et al. [ 27 ] and b) Chatterjee et al. [ 5 ] with predicted values from Eq. ( 13) using the 3 sets of parameters as given in Table 4. Solid lines are 
±20% margins. 


whole experimental data matrix, i.e. that which resulted from the flow 
of both sand and broken rice, then it was run with the sand and broken 
rice experimental data separately. The angle of repose and the (inverse 
of the) Hausner ratio are strongly dependent on the materials' proper¬ 
ties; for this reason, when determining the parameters with experi¬ 
mental data from either sand or broken rice, y and < were fixed to the 
values obtained from the whole experimental data matrix. Thus, three 
sets of parameters best fitting the experimental data using the present¬ 
ed model were obtained as given in Table 4. 


It should be noted that the aforementioned correlation for 
the MRT is applicable in the case of inclined rotary kilns, whether 
or not they are equipped with lifters or fitted with a dam at the 
outlet end. Using one of the proposed set of parameters in 
Table 4, the kiln must be operated with solid particles in cascading 
motion [25]. 

3.2.22. Predictive performance of the model. The performance of the 
model was assessed using the three sets of optimal parameters 


Chatterjee et al. [6] model, Eq.2 Shahhosseini et al. [4] model, Eq.5 



The authors model, Eq.13 



Fig. 9. Comparison of the experimental MRT with predicted values from Eqs. (2), (5) and (13) (using "sand & rice-fitted parameters”). 
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Table 5 

Characteristics of differing kilns with their operating parameters range of variation as given in the literature. 



L 

M 

D 

[m\ 

I/D 

H~ 

Materials 

Size 

0 

rT 

Pbulk 

Ptapped 

[rpm] 

S M 


Siift 

H 

[mm] 

[kg-m 3 ] 

[kg-m- 3 ] 

ri [kg-h- 1 ] 

[m 2 ] 

Bongo Njeng et al. [18] 

1.95 

0.10 

19.50 

Sand 

0.4-0.8 

39 

1422 

1543 

2-12 

2-5 0.69-2.5 

0-33.5 

0-0.008 





Rice 

1.8-4.4 

36 

889 

934 





Debacq [26] 

6.85 

0.75 

9.13 

U02F 2 

0.02-0.5 

42 

400 

500 

1.6-2.1 

1.43 508-895 

0 

0.438 





U 3 0 8 

0.01-0.05 

33 

1000 

1300 





Colin etal. [28] 

4.2 

0.21 

20.00 

Beech 

5-15 

42 

260 

284 

2-4 

1-2 4-8 

0 

0 

Saiet al. [27] 

5.9 

0.15 

39.33 

Umenite 

1.1 

27.4 

4200 

4410* 

1-3 

0.78-1.37 6-36 

15-28 

0 

Chatterjee et al. [5] 

2 

0.3 

6.67 

Iron ore 

3-6 

35 

2500 

2625* 

03-0.7 

1-3 15-60 

0-75 

0 


separately on the experimental results. For reasons of clarity, a criterion 
in performance assessments J was defined as follows: 



This criterion can be defined as the sum of the relative square errors 
between predicted and experimentally determined values of the MRT 
divided by the number of experiments. It will be used to determine 
which one among the three sets of optimal parameters gives the best 
predictions. The lower the value of this criterion, the better the 
predictions. 

The performance of the proposed model was first tested on the 
experimental data from this study as given in Fig. 5. Very good agree¬ 
ment was found between the experimental MRT obtained while pro¬ 
cessing sand and the predictions in general. However, as one may 
expect, the best predictions were obtained with the “sand-fitted pa¬ 
rameters” which displayed the lowest value of the criterion J as 
shown in Fig. 5a. Similar remarks can be made when looking at the 
experiments processing broken rice. There as well, good agreement 
was found between the experimental and predicted values of the 
MRT, and the “rice-fitted parameters" gave the least deviations 
from experimental values according to the criterion value (see 
Fig. 5b). In both cases the “sand & rice-parameters” give some good 
results, with values of the criterion sometimes very close to the min¬ 
imum observed. 

The predictive capacity of the proposed model was then tested on 
some experimental data from the literature. Kilns of different geome¬ 
tries and processing a range of materials at varying conditions were 
chosen. Only a few researchers investigating the residence time of par¬ 
ticles through rotary kilns have published their experimental data along 
with sufficient details on properties of the bulk materials. This is even 
truer for kilns equipped with lifters. The results from the work by 
Debacq [26] were used in this case. The results from kilns equipped 
with dams and no lifters were also tested. The latter are from the 
work by Sai et al. [27] and Chatterjee et al. [5], The results from the 
kiln used by Colin et al. [28], which was equipped with neither lifters 
nor exit dam, were also used in this study. The experimental results of 
the aforementioned kilns were used to evaluate the predictive perfor¬ 
mance of the model for kilns of varying geometries with differing 
types of material properties and operating conditions. Details about 
the geometry of the kilns and the properties of the materials are speci¬ 
fied in Appendix A. However, it should be pointed out that the mass 
flow rates achieved in most of these experiments extracted from the 
literature are (between 2 and 300 times) higher than those achieved 
in the present study. 

Fig. 6 shows a comparison of the experimental MRT reported in 
Debacq [26] with corresponding predictions. The model predictions 
are very good and the values of the criterion are very close for the 
three sets of parameters, but the “rice-fitted parameters” proved to be 


the best in this case with the lowest value of the criterion. The results 
are even better given that this industrial kiln was equipped with six 
rectangular lifters and processed a high mass flow rate of cohesive 
particles. 

The correlation was then used to calculate the MRT of beech chips 
flowing through the kiln used by Colin et al. [28], and operating without 
any internal fixtures or exit dam. It should be mentioned that this kiln is 
twice as big as the one used in the present study. Hence, these kilns have 
similar length-to-diameter ratios. The comparisons of predicted and ex¬ 
perimental values of the MRT are given in Fig. 7. Results show good 
agreement between the experimental and calculated MRT regardless 
of the set of parameters used. However, it should be noted that results 
obtained while using the “sand-fitted parameters" give the lowest 
value of the performance criterion followed by the “rice-fitted 
parameters.” 

Other tests were carried out with data from rotary kilns equipped 
with dams of differing heights at their outlet ends. These rotary kilns 
were those used by Sai et al. [27] and Chatterjee et al. [5] as previous¬ 
ly mentioned. Fig. 8 shows the results of the comparison of the ex¬ 
perimental and predicted values of the MRT. In the results obtained 
from the data of Sai et al. [27], the correlation gives a good predictive 
performance when using either the “rice-fitted parameters or sand & 
rice-fitted parameters," but revealed an even better agreement with 
the experimental MRT with the “sand-fitted parameters,” which 
gives a significantly low criterion compared with the other results. 
Chatterjee et al. [5] achieved a huge experimental matrix including 
5 types of dams, but their correlation for the MRT was not adapted 
for kilns fitted with exit dams (see Eq. (2)). The comparisons of the 
predicted and experimental values of the MRT show some devia¬ 
tions. In this case the predictive performance of the model seems to 
be somewhat diminished, as suggested by the elevated values ob¬ 
tained for criterion J. Nonetheless, about 30 to 50% of the results 
are within reasonable margins of the MRT experimentally measured. 
These poor results of the model can be attributed to (i) the difficul¬ 
ties in extrapolating to kilns with a rotational speed lower than 
2 rpm, leading most of time to small Froude numbers characteristic 
of beds in slipping motion [25], or perhaps (ii) the wide range ofvar- 
iation of the effective exit open diameter achieved by Chatterjee et al. 
[5] even if the results obtained in the case of Sai et al.'s [27] experi¬ 
ments were quite acceptable. 

To complete the analysis of the predictive capacity of the proposed 
model, other models developed for rotary kilns were applied in order 
to compare the MRT predictions. The correlation proposed by Chatterjee 
et al. [6] (see Eq. (2)) which does not account for the presence of lifters 
or an exit dam was used. The correlation of Shahhosseini et al. [4] (see 
Eq. (5)) which may apply to kilns equipped with lifters was also used. 
The predictions of these models are presented in Fig. 9, along with 
predictions obtained from Eq. (13) with the use of “sand & rice-fitted 
parameters.” The predictions from Chatterjee et al.'s [6] correlation un¬ 
derestimate the experimental MRT by almost half, while those from 
Shahhosseini et al.'s [4] correlation largely overestimate them. Further¬ 
more, the predictions from Eq. (5) display two patterns that are 
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associated with the two types of solid particles used in this study. These 
results give a large value of the performance criterion, whereas the pre¬ 
dictions from Eqs. (2) and (13) are more uniformly distributed, with a rel¬ 
atively lower value of the criterion. However, predictions from Eq. (2) 
display a vertical pattern linked to its incapacity to take into account the 
presence of lifters and exit dam. As shown in Fig. 9, predictions from 
Eq. (13) give the lowest value of the criterion compared to the other 
two models, failing to give good predictions of the experimental MRT. 

It appears, in accordance with the results of the present study, 
that the correlation proposed in this work for the MRT is also valid 
when applied to kilns with materials and operating conditions 
different from those used in the study. This confirms the excellent pre¬ 
dictive capacity of Eq. (13) compared to similar semi-empirical models 
which are often limited to a small range of conditions and operating 
variables. Regardless of the set of parameters, the predictions of the 
correlation are usually within reasonable margins of the experimental 
MRT. For that reason the authors recommend the use of the “sand & 
rice-fitted parameters.” 

4. Conclusions 

The effects of operating parameters, namely, the rotational speed, 
the kiln slope, the mass flow rate, the exit dam height, the type of mate¬ 
rials and shape of lifters, on the hold-up and mean residence time have 
been investigated. As far as these parameters are concerned, in this 
work it was observed that: 

• Increasing either the rotational speed or the kiln slope increased the 
kiln hold-up and the MRT of the charge. 

• Increasing the feed rate tended to increase the hold-up while the MRT 
of particles decreased. With an increase in the exit dam height the 
latter two increased slightly. 

• The use of rectangular lifters slightly increased the kiln hold-up 
and the MRT of particulate solids compared with straight or no 
lifters. 

• Sand and broken rice bulks have similar Hausner ratios though 
they have different bulk and tapped densities and different particle 
sizes. Their MRT were relatively close but the sand MRT were al¬ 
ways higher whereas the rice hold-ups were usually higher. 

A semi-empirical model was developed in order to predict the mean 
residence time of solid materials under steady state conditions in 
inclined rotating kilns, whether equipped or not with lifters and exit 
dams, as follows: 




The correlation consists of dimensionless factors accounting for 
(i) the basic operational parameters of the kiln, (ii) the properties of 
the bulk solids, (iii) the geometry of the kiln, (iv) the number, shape 
and presence of lifters, (v) the height of the exit dam. The parameters, 
namely k, a, j3, .... f have been determined from experimental data. 
The authors recommend the use of the “sand & rice-fitted parameters” 
as given in Table 4, which give good agreement with experimental 
data from this study as well as from published data. However, the sets 
of parameters proposed here are only applicable to beds in cascading 
(tumbling) motion as defined by Mellmann [25], Therefore, possible di¬ 
rections for further work could be (i) to extend the validity of the model 
to other types of bed motion with an adequate set of parameters, (ii) to 
enrich the experimental data matrix with results from the use of new 
types of solid materials to reaffirm the set of parameters determined. 


List of symbols 

D opl .„ Effective exit open diameter 

d p Average particle size 

D Kiln internal diameter 

RA&A Function 

F Solid feed rate 

factor Fitting parameters 

C Gas flow rate 

g Gravitational acceleration 

HU Hold-up 

k Fitting parameters 

J Performance criterion 

I Kiln length 

Ls Area occupied by solid particles contained 

in loading lifters 

M Mass flow rate 

MRT Mean residence time 

N Kiln rotational speed 

NL No lifters 

rinjt Total number of lifters 

Ns Total number of samples 

N t Total number of experiments 

RL Rectangular lifters 

rpm Rotation per minute 

RTD Residence time distribution 

S Kiln slope 

SL Straight lifters 

Sip Tube cross section not covered by solid 

particles in loading lifters 

Show i Section covered by solid particles in 

a lifter at horizontal position 
t Time 

t MRT 

Vf Volume of solids inside a lifter 

X hold-up volume fraction 
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Greek letters 


Fitting parameters 


At f Sampling times s 

8 Angle of repose 

p b ui k Bulk density kg-m~ * 2 3 4 

p p Particle density kg m 3 

ptme Tapped density kgm -3 

t Time of passage min 


Appendix A. Geometrical characteristics and operating conditions of 
some kilns 

Table 5 gives the characteristics of different kilns from the literature, 
along with the range of their operating conditions and characteristics of 
the materials used. For some materials the tapped density was missing 
in the published papers (those marked with * in Table 5). For these cases 
a literature review was done to evaluate the ability of the solid particles 
to flow. The flow ability is influenced by the particle size; fine particles 
(<0.1 mm) tend to be more cohesive and thus less free-flowing, where¬ 
as larger, denser particles tend to be free-flowing. It was decided to use 
a Hausner ratio of 1.05, characteristic of a free flowing bulk, for both 
Ilmenite and Iron ore. 
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